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Background.  Among  influenza  viruses,  type  A  viruses  exhibit  the  greatest  genetic  diversity,  infect  the  widest 
range  of  host  species,  and  cause  the  vast  majority  of  cases  of  severe  disease  in  humans,  including  cases  during  the 
great  pandemics.  The  hemagglutinin  1  (HA1)  domain  of  the  HA  protein  contains  the  highest  concentration  of 
epitopes  and,  correspondingly,  experiences  the  most  intense  positive  selection  pressure. 

Objectives.  We  sought  to  isolate  and  genetically  characterize  influenza  A  virus  subtype  H1N1  (A[H1N1]) 
circulating  in  Kenya  during  2007-2008,  using  the  HA1  protein. 

Methods.  Nasopharyngeal  swab  specimens  were  collected  from  patients  aged  >2  months  who  presented  to  8 
healthcare  facilities  in  Kenya  with  influenza-like  illness.  We  tested  specimens  for  seasonal  influenza  A  viruses, 
using  real-time  reverse-transcription  polymerase  chain  reaction  (RT-PCR).  Viruses  were  subtyped  using  subtype- 
specific  primers.  Specimens  positive  for  seasonal  A(H1N1)  were  inoculated  onto  Madin-Darby  canine  kidney  cells 
for  virus  isolation.  Viral  RNAs  were  extracted  from  isolates,  and  the  HA1  gene  was  amplified  by  RT-PCR,  fol¬ 
lowed  by  nucleotide  sequencing.  Nucleotide  sequences  were  assembled  using  BioEdit  and  translated  into  amino 
acid  codes,  using  DS  Gene,  version  1.5.  Multiple  sequence  alignments  were  performed  using  MUSCLE,  version 
3.6,  and  phylogenetic  analysis  was  performed  using  MrBayes  software. 

Results.  We  found  that,  similar  to  A/Brisbane/59/2007  (HlNl)-like  virus,  which  was  included  in  the  southern 
hemisphere  vaccine  for  the  2009  influenza  season,  all  2007  Kenyan  viruses  had  D39N,  R77K,  T132V,  K149R,  and 
E277K  amino  acid  substitutions,  compared  with  A/Solomon  Islands/3/2006  (HlNl)-like  virus,  a  component  of 
the  southern  hemisphere  vaccine  for  the  2008  influenza  season.  However,  the  majority  of  2008  viruses  from 
Kenya  also  had  R192K  and  R226Q  substitutions,  compared  with  A/Solomon  Islands/3/2006  (HlNl)-like  virus. 
These  2  changes  occurred  at  the  receptor  binding  site.  The  majority  of  the  2008  Kenyan  isolates  contained  N187S, 
G189N,  and  A193T  mutations,  which  differed  from  A/Brisbane/59/2007  (HlNl)-like  virus.  The  A193T  substitu¬ 
tion  is  involved  in  binding  the  sialic  acid  receptor.  Phylogenetically,  the  2008  Kenyan  isolates  grouped  into  2 
clusters.  The  main  cluster  contained  viruses  with  N187S  and  A193T  changes;  residue  187  is  involved  in  receptor 
binding,  whereas  residue  193  is  at  antigenic  site  Sb. 

Conclusion.  Overall,  the  major  genetic  variations  that  occurred  in  seasonal  A(H1)  viruses  either  affected  re¬ 
ceptor  binding  or  altered  epitopes  at  the  immunodominant  sites.  These  genetic  variations  in  seasonal  A(H1N1) 
isolated  in  Kenya  during  2007-2008  highlight  the  importance  of  continuing  surveillance  and  characterization  of 
emerging  drift  variants  of  influenza  virus  in  Africa. 
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Among  influenza  viruses,  type  A  viruses  cause  the  vast  major¬ 
ity  of  severe  disease  in  humans,  including  pandemics.  They 
also  infect  the  widest  range  of  host  species  and  exhibit  the  great¬ 
est  genetic  diversity  [1],  The  genetic  diversity  of  influenza 
viruses  is  caused  by  the  error-prone  RNA  polymerase  during 
genome  replication  [2,  3].  The  error-prone  replication  and  short 
generation  times  result  in  large,  highly  diverse  replicating  popu¬ 
lations  of  viral  quasi-species  upon  which  natural  selection  oper¬ 
ates  [4].  Investigating  the  complex  factors  that  affect  this  natural 
selection  is  key  to  understanding  influenza  A  virus  evolution. 
Modern  molecular  and  phylogenetic  analysis  tools  are  often 
applied  to  understand  these  processes  [5, 6], 

While  both  hemagglutinin  (HA)  and  neuraminidase  (NA) 
proteins  contain  immunogenic  sites  at  which  immune-driven 
natural  selection  happens,  the  N-terminal  (HA1)  domain 
of  the  HA  protein  contains  the  maximum  concentration  of 
epitopes  and  experiences  the  strongest  positive  selection  pres¬ 
sure  [7],  The  majority  of  neutralizing  antibodies  are  directed 
toward  5  distinct  immunogenic  sites  on  the  HA1  globular 
domain:  Sa,  Sb,  Cal,  Ca2,  and  Cb  [8,  9].  Each  of  these  immu¬ 
nogenic  sites  is  a  conformational  entity  formed  by  a  number 
of  amino  acid  residues  folded  in  3-dimensional  space  [8,  9], 
Any  change  occurring  at  these  immunogenic  sites  will  affect 
the  ability  of  preexisting  antibodies  to  neutralize  the  virus. 
Because  the  human  immune  response  to  viral  infection  is 
incompletely  cross-protective,  amino  acid  variants  of  the 
HA1  protein  that  allow  the  virus  to  evade  immune  response 
to  infect  more  hosts  are  favored  by  natural  selection  pro¬ 
cesses  [7].  In  addition  to  the  central  role  of  the  HA1  domain 
as  an  immune  response  determinant,  it  is  also  the  influenza 
virus  ligand  containing  the  binding  site  for  attachment  to  the 
host  cell  receptors  [10,  11].  The  high  affinity  of  HA1  for 
human-like  a(2,6)  receptors  and  the  significant  difference  in 
sequence  and  immunogenicity  of  HA1  among  existing  season¬ 
al  and  vaccine  strains  are  the  hallmark  pathogenicity  determi¬ 
nants  for  influenza  viruses,  influencing  the  ease  of  human- 
to-human  transmission.  Within  HA1,  D190,  and  D225  residues 
are  key  determinants  for  effective  binding  to  human-like  a(2,6) 
receptors  [10,  12,  13].  A  single  mutation,  D225G,  reduces  HA’s 
binding  affinity  for  a(2,6)  receptors  [12,  13]  and  the  infectivity 
of  the  virus  [10],  while  a  double  variant,  D190E/D225G,  elimi¬ 
nates  HA’s  binding  affinity  to  a(2,6)  receptors  [12,  13]  and 
renders  the  virus  noninfectious  [10]. 

Kenya  is  an  East  African  country  that  straddles  the  equator 
and  has  a  tropical  climate.  Although  2  reports  showing  that 
influenza  A  virus  subtype  H3N2  and  influenza  B  virus  circu¬ 
lated  in  Kenya  between  2004  and  2006  [14,  15],  there  has  not 
been  documentation  of  the  presence  and  molecular  characteri¬ 
zation  of  influenza  A  virus  subtype  H1N1  (A[H1N1])  in 
Kenya.  This  is  mainly  due  to  a  lack  of  influenza  surveillance 
in  the  country  [14,  15].  Although  influenza  vaccine  is  rarely 
used  in  Kenya,  the  Kenya  Ministry  of  Health  recommends  the 


use  of  the  southern  hemisphere  influenza  vaccine  formulation 
for  its  citizens  and  other  residents. 

In  2006,  in  collaboration  with  the  US  Centers  for  Disease 
Control  and  Prevention  (CDC)  and  the  Kenya  Ministry  of 
Public  Health  and  Sanitation,  the  US  Army  Medical  Research 
Unit  established  syndromic  and  virologic  surveillance  for  influ¬ 
enza  at  sentinel  health  facilities  throughout  Kenya  [14].  From 
respiratory  specimens  collected  through  this  surveillance  system 
during  2007-2008,  we  isolated  and  analyzed  A(H1N1).  By  using 
the  nucleotide  sequences  of  the  globular  HA1  region  of  the  HA 
gene,  we  provide  the  first  description  of  the  molecular  and  phy¬ 
logenetic  characteristics  of  A(H1N1)  circulating  in  Kenya. 

MATERIALS  AND  METHODS 

Setting 

We  conducted  surveillance  for  patients  with  influenza-like 
illness  (ILI)  at  1  provincial  hospital  in  Kisumu  city  in  western 
Kenya  and  in  7  district  hospitals  in  the  cities  of  Nairobi, 
Mombasa,  Malindi,  Isiolo,  Kericho,  Kisii,  and  Alupe.  Surveil¬ 
lance  hospitals  were  chosen  because  they  were  geographically 
and  demographically  representative  of  the  population  of  Kenya. 

Study  Population 

Patients  with  ILI  included  anyone  >2  months  of  age  present¬ 
ing  to  a  surveillance  hospital  outpatient  clinic  with  fever  and 
either  cough  or  sore  throat.  We  excluded  presentations  of  ILI 
in  patients  with  exudative  pharyngitis  or  tonsillitis  or  symptom 
onset  >72  hours  prior  to  presentation.  Surveillance  officers  en¬ 
rolled  the  first  5  case  patients  with  ILI  presenting  to  each  sur¬ 
veillance  hospital  daily  and  obtained  informed  consent  from 
patients  for  study  participation.  The  Kenya  Medical  Research 
Institute  (KEMRI),  the  Walter  Reed  Army  Institute  of  Research 
(WRAIR),  and  the  CDC  institutional  review  boards  reviewed 
and  approved  the  study  protocol  under  protocol  numbers  SSC 
981,  WRAIR  1267,  and  CDC  4716,  respectively. 

Specimen  Collection 

Surveillance  officers  collected  duplicate  nasopharyngeal 
samples  from  each  case  patient,  using  Dacron-flocked  swabs. 
Swabs  were  placed  in  a  1-mL  cryovial  containing  virus  trans¬ 
port  medium,  kept  at  4°C,  and  stored  in  a  liquid  nitrogen  dry 
shipper  within  8  hours  of  collection.  All  samples  were  trans¬ 
ported  from  the  surveillance  sites  to  the  Kenya  National  Influ¬ 
enza  Center  (NIC)  laboratory  within  1  week.  The  cold  chain 
was  maintained  throughout. 

Laboratory  Testing 

Real-Time  Reverse-Transcription  Polymerase  Chain 
Reaction  (RT-PCR) 

For  detection  of  influenza  virus  in  patient  samples,  total  RNA 
was  extracted  from  the  samples,  using  the  QIAamp  Viral  RNA 
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Mini  Kit  (Qiagen,  Valencia,  CA)  according  to  the  manufactur¬ 
er’s  instructions.  One-step  real-time  PCR  was  performed  on 
an  ABI  7500  Fast  platform,  according  to  the  CDC  protocol  for 
the  identification  of  A(H1)  viruses.  The  total  reaction  volumes 
were  25  pL,  containing  0.5  pL  of  superscript  Ill/Platinum 
enzyme  mix  (Invitrogen,  Paisley,  United  Kingdom),  5.5  pL  of 
H20,  12.5  pL  of  2X  buffer,  0.8  pM  forward  primer  (40  pM), 
0.8  pM  reverse  primer  (40  pM),  and  0.2  pM  probe  (10  pM). 
Each  RNA  sample  was  tested  for  6  sets  of  gene  targets:  matrix 
gene  segment,  to  identify  influenza  A  and  B  viruses;  HA  gene 
segment,  to  subtype  HI,  H3,  and  H5  viruses;  and  human  RNP 
gene  segment,  to  test  the  RNA  extraction  procedure.  Reverse 
transcription  was  achieved  at  50°C  for  30  minutes  and  95°C. 
PCR  was  achieved  after  45  cycles  of  95°C  for  15  seconds  and 
55°C  for  30  seconds. 

Virus  Isolation 

Influenza  virus  isolation  was  performed  at  the  NIC  in  Madin- 
Darby  canine  kidney  cells.  Afterward,  we  conducted  a  hemag¬ 
glutination  inhibition  (HAI)  assay,  using  guinea  pig  red  blood 
cells  and  reference  antiserum  in  accordance  with  CDC  proto¬ 
cols.  Inoculated  cells  were  incubated  at  37°C  with  5%  C02  and 
observed  daily  for  10  days  for  visual  cytopathic  effect,  using 
an  inverted  microscope.  When  the  cytopathic  effect  was  ob¬ 
served,  the  supernatant  fluid  was  collected,  and  the  hemagglu¬ 
tination  titer  was  measured.  HAI  testing  was  performed  on 
high-titer  (>32:1)  samples.  Isolates  with  a  low  hemagglutina¬ 
tion  titer  (<32:1 )  were  repassaged  once.  On  the  basis  of  the 
variability  in  HAI  titers,  we  selected  46  isolates  for  further 
molecular  genotyping,  and  these  were  confirmed  as  A(H1N1). 
The  guinea  pigs  used  to  provide  blood  for  HA  and  HAI  assays 
of  this  report  were  reared  in  accordance  with  KEMRI  guide¬ 
lines  on  animal  care. 

PCR  Assay  and  Nucleotide  Sequencing  of  the  HAI 
Gene  Segment 

A(H1N1)  isolates  were  sent  to  US  Air  Force  School  of  Aero¬ 
space  Medicine  (USAFSAM;  San  Antonio,  TX),  where  they 
were  processed  for  genome  sequencing.  The  USAFSAM 
served  as  the  confirmatory  laboratory  for  HAI  tests  done  at 
KEMRI.  Isolates  were  molecularly  typed  and  subtyped  at  the 
USAFSAM,  and  results  were  forwarded  to  the  CDC  in 
Atlanta,  a  WHO  Influenza  Collaborating  Centre. 

RNA  was  extracted  from  48-hour  shell  vial  cultures  at  the 
USAFSAM,  using  the  M48  BioRobot  and  Mag  Attract  Virus 
Mini  M48  Kit  (Qiagen)  according  to  the  manufacturer’s  pro¬ 
tocols.  For  RT-PCR  amplification,  5  pL  of  RNA  was  added  to 
a  50-pL  master  mixture  containing  2  x  reaction  buffer, 
1  x  enzyme  mixture,  and  400  nmol/L  each  of  CDC  primers 
(H1HA_26_F_M13,  5'-TGT  AAA  ACG  ACG  GCC  AGT  CAA 
CCA  AAA  TGA  AAG-3';  and  H1HA_1238_R_M13,  5'-CAG 


GAA  ACA  GCT  ATG  ACC  AGC  TGT  GAA  TTG  RGT  GTT 
CAT  TT-3'),  using  the  Superscript  III  One-Step  RT-PCR 
System  (Invitrogen,  Carlsbad,  CA).  RT-PCR  thermocycling 
consisted  of  cycles  at  50°C  for  30  minutes,  hot  start  activation 
at  94°C  for  2  minutes,  40  amplification  cycles  at  94°C  for  15 
seconds,  at  52°C  for  30  seconds,  and  at  68°C  for  75  seconds, 
with  a  final  extension  cycle  at  68°C  for  5  minutes.  All  PCR 
products  were  visualized  as  1200-base  pair  bands  after  gel 
electrophoresis  in  2%  precast  agarose  gels  stained  with  ethid- 
ium  bromide  (Invitrogen)  under  UV  illumination.  PCR  prod¬ 
ucts  were  purified  using  the  QIAquick  PCR  Purification  Kit 
(Qiagen).  Incorporation  of  fluorescent  nucleotides  was  facili¬ 
tated  using  the  Big  Dye  Terminator  v3.1  Kit  (Applied  Biosys¬ 
tems,  Foster  City,  CA)  and  sequenced  using  terminal 
H 1 H  A_26_F_M  1 3  and  H1HA_1238_R_M13  RT-PCR 
primers  and  6  internal  oligonucleotides:  H1HA_522_R_M13 
(5'-CAG  GAA  ACA  GCT  ATG  ACC  CGT  CAG  CCA  TAR 
CAA  ATT  T-3'),  H1HA_239_F_M13  (5'-TGT  AAA  ACG 
ACG  GCC  AGT  TGG  GTA  AYT  GCA  GCR  TTG  C-3'), 
H1HA_760_R_M13  (5'-CAG  GAA  ACA  GCT  ATG  ACC 
GTC  CAG  TAG  TAR  TTR  ATT  CT-3'),  H1HA_F3_487  (5'- 
TGT  AAA  ACG  ACG  GCC  AGT  AAT  TTG  CYA  TGG  CTG 
ACG-3'),  H1HA_1063_R_M13  (5'-CAG  GAA  ACA  GCT 
ATG  ACC  GCA  ATG  GCT  CCA  AAC  AAA  CCT  CT-3'),  and 
H1HA_F4_682  (5'-TGT  AAA  ACG  ACG  GCC  AGT  CCA 
GAA  ATA  RCC  AAA  AGA  C-3'). 

HAI  Nucleotide  Genetic  Analyses 

Nucleotide  contigs  were  assembled  using  the  Contig  Assembly 
program  in  Bioedit  [16].  Protein  translations  were  performed 
using  DS  Gene,  version  1.5  (Accelrys,  Cambridge,  United 
Kingdom).  Multiple  sequence  alignments  were  performed  with 
MUSCLE,  version  3.6  [17].  Phylogenetic  analyses  were  per¬ 
formed  with  MrBayes  3.1  [18].  The  generalized  time-reversible 
model  parameters  and  priors  were  incorporated  into  the  nexus 
file  for  execution  in  MrBayes.  The  data  were  executed  in 
MrBayes  by  running  1  million  Monte  Carlo  Markov  chains 
with  a  sampling  frequency  of  1000.  The  phylogenetic  tree  was 
visualized  using  FigTree,  version  1.3.1  [19].  The  protein  se¬ 
quences  for  the  Kenyan  A(H1N1)  isolates  are  available  from 
GenBank  (accession  numbers  FJ662910.1-FJ662927.1, 
FJ265654.1-FJ265662.1,  and  EU204969.1-EU770350.1). 

The  nucleotide  and  deduced  amino  acid  sequences  of  HAI 
from  isolates  were  compared  with  those  of  vaccine  strains,  as 
well  as  with  those  of  strains  from  another  African  country 
(South  Africa)  and  2  other  tropical  countries  (Malaysia  and 
Thailand).  Thus,  the  HAI  of  A/Brisbane/59/2007  (HlNl)-like 
virus  (the  recommended  southern  hemisphere  vaccine  strain 
for  the  2009  influenza  season),  A/Solomon  Islands/3/2006 
(HlNl)-like  virus  (the  recommended  southern  hemisphere 
vaccine  strain  for  the  2008  season),  and  A/New  Caledonia/20/ 
99(H1N1)-Iike  virus  (the  recommended  southern  hemisphere 
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vaccine  strain  for  the  2000-2007  seasons)  were  included  in  the 
analysis. 

RESULTS 

From  January  2007  through  November  2008,  we  enrolled  and 
collected  nasopharyngeal  swab  specimens  from  5898  children 
and  adults  who  received  a  diagnosis  of  ILI.  The  median  age  of 
enrolled  case  patients  was  18  months  (range,  2  months-75 
years),  and  5629  specimens  (94.5%)  were  from  children.  None 
of  the  patients  sampled  had  previously  been  vaccinated 
against  influenza. 

All  samples  tested  by  real-time  RT-PCR  were  positive  for 
the  human  RNP  housekeeping  gene,  but  only  1177  samples 
were  positive  for  the  matrix  gene  of  influenza  A  virus.  Of  the 
1177  matrix  gene-positive  samples,  72  representative  isolates 
(6.1%)  for  each  month  of  collection  and  from  diverse  geo¬ 
graphical  regions  of  the  country  were  randomly  selected  for 
further  testing.  On  the  basis  of  variability  in  HAI  titers,  46  of 
72  isolates  (65%)  were  confirmed  as  A(H1N1)  and  were  select¬ 
ed  for  further  molecular  genotyping. 

Compared  with  the  reference  influenza  vaccine  strains  and 
the  influenza  strains  from  South  Africa  and  tropical  countries 
(Malaysia  and  Thailand),  all  2007  Kenyan  A(H1N1)  isolates 
were  most  similar  (98.6%-100%,  based  on  nucleotides; 
and  98.1%-100%,  based  on  amino  acids)  to  A/Brisbane/59/ 
2007  (HlNl)-like  virus,  followed  by  A/Solomon  Islands/3/ 
2006  (HlNl)-like  virus  (96.7%-99.9%,  based  on  nucleo¬ 
tides;  and  96.5%-100%,  based  on  amino  acids),  and  A/New 


Caledonia/20/99(H1N1)-Iike  virus  (95.9%-99.9%,  based  on 
nucleotides;  and  98.1%-99.7%,  based  on  amino  acids). 

Mutational  analyses  of  the  HAI  protein  that  show  evolution 
of  the  2007-2008  Kenyan  isolates  is  presented  in  Table  1. 
Compared  with  A/Solomon  Island/3/2006  (HlNl)-like  virus, 
all  2007-2008  Kenyan  isolates  had  D39N,  R77K,  T132V, 
K149R,  and  E277K  substitutions.  In  addition,  42%  of  the  2007 
isolates  had  the  Q108H  change,  while  31.5%  of  the  isolates 
from  the  same  year  had  the  S125N  mutation.  We  observed 
that  2  of  the  2007  isolates  had  D225G  and  that  all  of  the  2008 
isolates  had  R226Q  changes.  Furthermore,  of  the  2008  isolates, 
6  (21%),  26  (93%),  and  8  (29%)  had  N187S,  R192,K  and 
A193T  changes  respectively,  and  24  (86%)  of  the  2008  Kenyan 
isolates  had  a  change  affecting  residue  189.  Of  these,  23  (96%) 
were  G189N,  and  1  isolate  had  a  G189A  mutation. 

All  of  the  2007  Kenyan  isolates  had  the  Q108H,  S125N/ 
S125R,  and  D225G  changes  in  exactly  the  same  proportions  as 
A/Solomon  Islands/3/2006  (HlNl)-like  virus.  Q108H  and 
S125N  are  parallel  mutations  found  in  Kenyan  isolates,  but 
have  not  been  seen  anywhere  else  globally,  during  and  after 
the  study  period.  In  addition,  7  (37%)  of  the  2007  Kenyan  iso¬ 
lates  had  a  mutation  at  residue  190.  The  majority  (6  [86%])  of 
these  had  a  D190N  change,  and  only  a  single  isolate  had  the 
previously  reported  D190G  mutation.  The  majority  (24 
[86%])  of  the  2008  Kenyan  viruses  had  a  mutation  at  residue 
189;  23  (96%)  were  G189N  and  the  remaining  had  a  G189A 
mutation.  Furthermore,  we  discovered  that  of  the  2008 
Kenyan  isolates,  6  (21%)  and  8  (29%)  had  N187S  and  A193T 
mutations,  respectively. 


Table  1.  Major  Amino  Acid  Changes  Present  within  the  HAI  Protein  in  2007-2008  Kenyan  Influenza  A  (H1N1)  Isolates  Compared  to  A/ 
Solomon  lslands/3/2006  and  A/Brisbane/59/2007  Vaccine  Strains 

Year 

Vaccine  Strain 

2007 

2008 

A/Solomon  lslands/3/2006 

D39N  (19/19) 

D39N  (28/28) 

R77K  (1 9/1 9) 

R77K  (28/28) 

Q108H  (8/19) 

T132V  (28/28) 

S125N  (6/19);  SI  21 R  (1/1 9) 

K149R  (28/28) 

T1 32V  (19/19) 

N187S  (6/28)  [Sb] 

K149R  (19/19) 

G189N  (23/28);  G189A  (1/28)  [Sb] 

D225G  (2/19)  [Ca2;  SA1 

R192K  (26/28)  [SA,  Sb] 

E277K  (19/1 9) 

A193T  (8/28)  [SA,  Sb] 

R226Q  (28/28)  [SA] 

E277K  (28/28) 

A/Brisbane/59/2007 

Q108H  (8/19) 

N187S  (6/28)  [Sb] 

S125N  (6/19);  S125R  (1/19) 

G189N  (23/28);  G189A  (1/28)  [Sb] 

D190N  (6/19);  D190G  (1/19)  [Sb,  SA] 

K192R  (2/28)  [Sb,  SA] 

D225G  (2/19)  [Ca2,  SA] 

A193T  (8/28)  [SA,  Sb] 

The  numbers  in  parentheses  are  proportions  of  the  isolates  with  the  indicated  change.  Ca2,  Sb,  and  SA  denote  occurrence  at  antigenic  sites  Ca2  and  Sb  and 
binding  to  sialic  acid  receptor,  respectively. 
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A_New  Caledonia_20_99  H1N1 

- A_Madagascar_57794_2000_H1N1 

- [M 


- A_Johannesburg_28_05_H1N1 

A_Johanr»esburg_ 1 46  05  H 1 N 1 


100% 


J39N.  R77K.  T132V,  K149R  E277> 


R188K 


•  A_Sok>monJstands_3_2006_H1N1 


A  lsiolo_5880_2008  H1N1 
A  Mbagathi  484 1  2007  H1N1 
A  Mba 


■  A_Malaysia_330_2007_H1N1 


- A  Mbagathi  < 

•  A  Kisumu_4830_20O7  HINT 


ithi  4833_2007_H INI 


-  AKericho  4822  2007  H1N1 
■  ATtfalindi  4836_2007_RlN1 

— - - A_ Johannesburg  21  2008_H1N1 

-  A  Kisii  5890  2608  FlINI 

•  A  Malindi  4838  20(57  H1N1 


-  A_Kencho  5883  2008  RlNI  ' 

-  A_Kisti_4824  2007  HfNl 

-  A  Mbagathi  4819  2007  H1N1 

-  A3nsbane39 2d07  HfNl 


-  A  Thailand_ 39_2008_H  1 N 1 


-  A_Thailand_530_2008_H  1 N 1 

•  A  Johannesburg  67  2007  H1N1 
•  A'  Thailand  309  2008  HfNl 
•  A^Kisii  4823  2007  H1N1 


-  A  Malindi_4837l200r  HINl 
■  AjJohannesburg  28  2008  H1N1 

-  A_Cape_Tovvn  28_20O8  RlNI 
- A  Thailand  50f  2007  H1N1 

-  A'Thailand  500  2007  H1N1 


HE 


—  A  Mbagathi  4640  2007  R1N1 
■  A  MEagalhi  4826  2007  HT1M1 

,  Kericho  4831  "2007  H1N1 


67% 


•  AMbagathT  4827  2007  H1N1 
-  A  Mbagathf  4844"  200 T  H1N1 

•  A  "Kisii  4839  2007  H1M1 
■  A2Malindi_4835_2007_H1N1 

-  A  Mbagathi_4825_2007  H1N1 

-  A“Mbagalhi  4834  2007 “MINI 
-“A  Kencho  7535  2008  H1N1 


•  A  Kencho  7525  2008  H1N1 


-  A  Kisii  7565  2008  “H1N1 

-  A  "Kisumu  5933  2008  H1N1 
■  A  Kisumu  '6543  "2008  >I1N1 

■  A  Kenya  6901  2008fH1N1 


A  Kisumu  761 1  2008  HTlNl 
A  Kisii  7554  2008  HfNl 
A  Kisif  7543"  2008  H1N1 
A _KehcFio_7530_2008_H  1 N 1 
A_Kericho_7529l2008_H  1 N 1 
A  Kisumu  5923  2008  MINI 
A  Kisii  5696  2008  HTN1 
A_Kisii_5886_2008_H1N1 

hntin - A  Kisumu  5928  2008  H1N1 

-  A  Kisumu  5934  2006_H1N1 

A  Mbagathi  7596_2008_H1N1 
A  Kisii  7577  2008  “H INI 
”  "  A_Th«Mn3_778_20O8  MINI 

A  Johannesburg_35  20O8  H1N1 
A  Johannesburg'  s/* 2008  H1N1 
A  Kisii  7576  2008  HINf 
A  Kisli  “7570  "2008  “Hi  N1 
A  Kisii  7562  2008  H1N1 
A  Kisii “7559  2008  H1N1 
A  Kisif  7547“ 2008  H1N1 
A  T<isiC754 1J2008  “Hi  N 1 
A_Johannesburg_25  2008  H1N1 
A  Johannesburg “1 0 “2008 “HI  N 1 
A  Johannnafcurg  2/9  2008  H 1 N 1 

rair - A  Johannesburg  46  2008  MINI 

—  A  Johannesburg  34  2008  H1N1 


00080 


Figure  1.  Phylogenetic  relationships  of  HA1  proteins.  The  amino  acid  chains  for  the  HA1  region  of  2007-2008  Kenyan  isolates  and  vaccine  strains  for  2000-2009  and  representative  South  African  and 
tropical  (Malaysian  and  Thai)  strains  of  this  period  were  analyzed.  Online,  green  represents  vaccines  strains,  blue  represents  the  2007  isolates,  and  red  represents  the  2008  isolates.  Black  taxa  are  earlier 
strains  from  South  Africa.  Labels  on  the  branches  represent  the  main  amino  acid  changes  for  fixed  mutations  defining  the  main  clades.  The  tree  is  unrooted,  and  all  horizontal  branch  lengths  are  drawn  to  a 
scale  of  substitutions  per  site  (as  shown  by  the  scale  bar).  The  percentages  represent  posterior  probabilities  of  each  branch  representing  clade  credibility. 
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A  phylogenetic  comparison  of  the  2007-2008  Kenyan  iso¬ 
lates  to  the  2008  and  2009  vaccine  strains  and  the  contempora¬ 
neously  circulating  South  African,  Malaysian,  and  Thai  isolates 
is  presented  in  Figure  1.  All  2007  Kenyan  isolates  clustered  in 
a  single  clade  together  with  A/Brisbane/59/2007  (HlNl)-like 
virus  but  away  from  A/Solomon  Islands/3/2006  (HlNl)-like 
virus.  Within  this  clade,  many  of  the  Kenyan  viruses  formed  a 
subbranch  away  from  A/Brisbane/59/2007  (HlNl)-like  virus, 
characterized  by  Q108H  and  S125N  changes.  A  majority 
(82%)  of  the  2008  Kenyan  isolates  had  a  G189N  mutation, 
which  was  not  present  in  A/Brisbane/59/2007  (HlNl)-like 
virus,  and  formed  a  clade  distinct  from  the  2007  Kenyan  iso¬ 
lates.  This  clade  bifurcated  further  into  2  branches:  the  first 
contained  viruses  without  the  A193T  change,  and  the  second 
contained  viruses  that  had  this  mutation.  The  majority  of  the 
viruses  present  in  the  A193T  branch  had  the  N187S  mutation. 
All  Kenyan  viruses  on  the  N187S  branch  were  derived  from  a 
single  sentinel  site,  and  only  viruses  isolated  from  Johannes¬ 
burg  in  2008  were  also  present  on  this  branch. 

DISCUSSION 

We  have  documented  the  evolution  of  A(H1N1)  viruses  in 
Kenya  during  a  23-month  period  in  2007-2008.  Compared 
with  the  2  vaccine  strains,  A/Solomon  Island/3/2006  (H1N1)- 
like  virus  and  A/Brisbane/59/2007  (HlNl)-like  virus,  the 
2007-2008  Kenyan  viruses  had  HA1  protein  mutations  at  5 
important  ligand  interacting  positions:  190,  192,  193,  225,  and 
226.  These  residues  are  involved  in  binding  the  human-like 
a(2,6)  ligand  [10,  12,  13].  The  mutations  detected  include 
D190N/G,  R192K  <->K192R,  A193T,  D225G,  and  R226Q. 
Amino  acids  at  position  190,  192,  and  193  contact  the  exten¬ 
sion  region  of  a(2,6)  ligand,  and  the  residue  at  226  interacts 
with  the  base  region  of  the  ligand  [12].  A  single  D225G  muta¬ 
tion  reduces  binding  affinity  for  a(2,6)  receptors  [12,  13]  and 
lowers  the  infectivity  of  the  virus  [10].  Double  D190E/D225G 
mutations  are  lethal  for  the  virus  as  they  abrogate  binding  to 
a(2,6)  receptors  [12,  13],  making  the  virus  noninfectious  [10], 
Furthermore,  mutations  occurring  at  positions  190,  225,  and 
226  were  present  at  the  Ca2  antigenic  site. 

Overall,  the  ligand-interacting  mutations  we  observed  in 
Kenyan  viruses  are  cumulatively  deleterious  to  the  virus  and 
individually  represent  viruses  evolving  to  less  fit  variants.  Con¬ 
sidering  the  relatively  large  numbers  of  such  viruses  we  ob¬ 
served  among  the  Kenyan  isolates,  A(H1N1)  viruses  seem  to 
be  evolving  at  high  rates  in  Kenya,  leading  to  this  dispropor¬ 
tionate  level  of  quasi-species  among  the  isolates. 

We  further  sought  to  investigate  in  silico  the  effects  that  mu¬ 
tations  occurring  in  the  HA1  protein  would  have  on  the  im- 
munogenicity  of  the  viruses.  We  found  5  mutations — D39N, 
R77K,  T132V,  K149R,  and  E277K — that  were  common  to  all 
viruses  circulating  in  Kenya  during  the  study  period.  None  of 


these  mutations  was  present  at  an  antigenic  or  receptor 
binding  site.  Similarly  mutations  at  positions  192  and  193  were 
present  at  the  Sb  antigenic  site  [8,  9].  In  addition,  we  observed 
that  the  majority  of  the  2008  Kenyan  viruses  had  changes  at 
positions  187  and  189,  both  located  at  antigenic  site  Sb.  The 
G189N  change  found  in  82%  of  the  2008  viruses  is  a  major 
change  involving  replacement  of  the  tiny  nonpolar  glycine 
residue  with  the  charged  bulkier  asparagine.  Because  of  steric 
and  electrostatic  interactions,  this  change  is  expected  to  alter 
this  epitope  to  abrogate  binding  of  any  preexisting  antibodies 
to  the  Sb  site.  In  general,  the  most  commonly  mutated  antigen¬ 
ic  site  among  the  2007-2008  Kenyan  viruses  was  the  Sb  site. 
The  apparent  positive  selection  on  the  Sb  antigenic  site  is  con¬ 
sistent  with  previous  reports  that  have  shown  that  this  site  is 
under  much  higher  pressure  for  mutations  than  the  others  [9] . 

All  2008  Kenyan  isolates  had  changes  at  the  receptor 
binding  site  and  multiple  mutations  affecting  antigenic  sites 
and  were  genetically  distinct  as  compared  to  A/Solomon 
Island/3/2006  (HlNl)-like  virus.  Similar  changes  were  seen 
globally  and  led  to  the  WHO  to  recommend  changing  the 
vaccine  A(H1)  strain  for  2009  to  A  Brisbane/59/2007 
(HlNl)-like  virus.  Overall,  our  data  show  that  the  majority  of 
the  2007  Kenyan  isolates  did  not  evolve  to  become  genetically 
different  from  the  vaccine  strain  and  that  the  southern  hemi¬ 
sphere  vaccine  was  therefore  well  matched  to  viruses  circulat¬ 
ing  in  Kenya  in  2007. 

The  clustering  on  the  same  phylogenetic  branch  as  A/Bris¬ 
bane/59/2007  (HlNl)-like  virus  for  all  the  Kenyan  2007  iso¬ 
lates,  in  contrast  to  the  single  Malaysian  2007  isolate,  which 
clustered  on  the  same  branch  as  A/Solomon  Islands/3/2006 
(HlNl)-like  virus,  suggested  that  the  Kenyan  2007  isolates 
were  closely  related  to  viruses  circulating  in  other  parts  of  the 
world  in  2007.  On  the  other  hand,  clustering  of  the  2008 
Kenyan  isolates  into  2  branches  characterized  by  N187S  and 
A193T  changes  and  on  a  separate  branch  from  A/Brisbane/ 
59/2007  (HlNl)-like  virus  showed  that  these  viruses  were 
genetically  variant  from  the  A(H1)  component  of  the  2009 
southern  hemisphere  vaccine.  These  observations  indicate  that 
the  viruses  that  circulated  in  Kenya  in  2008  drifted  away  from 
the  vaccine  strain  and  that  the  different  clusters  could  be  due 
to  different  introduction  of  viruses  in  the  country. 

Although  genetic  analysis  of  the  HA1  gene  is  important  in 
informing  our  understanding  of  influenza  virus  evolution  and 
antigenic  changes,  the  lack  of  hemagglutination  inhibition 
(HAI)  data  for  the  Kenyan  viruses  is  a  major  limitation  in  this 
study.  Furthermore,  although  HA  is  the  principal  antigen  that 
elicits  host  immune  response  and  is  responsible  for  binding  to 
sialic  acid  receptors  thus  mediating  viral  entry  into  host  cells 
[20],  HA  is  not  the  only  protein  against  which  antibodies  are 
directed.  Other  proteins,  including  NA,  play  a  role  in  overall 
host  humoral  response  to  the  virus.  However,  antibodies  to 
NA  do  not  neutralize  the  virus  as  do  HA  antibodies  [21]. 
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Thus,  analyses  of  HA1  genomic  data  alone  are  routinely 
applied  to  provide  preliminary  understanding  of  influenza 
viruses  [22-24].  Another  limitation  of  this  study  is  the  se¬ 
quencing  of  isolates  instead  of  viruses  from  clinical  specimens. 
This  may  have  introduced  mutations  not  found  in  the  original 
clinical  samples. 

Our  study  is  the  first  to  describe  the  genetic  variations  of  A 
(H1N1)  found  in  Kenya.  The  findings  regarding  changes  at  the 
receptor  binding  sites  suggest  that  in  Kenya,  A(H1N1)  is  under¬ 
going  directional  evolution  in  the  HA  protein  as  has  been  ob¬ 
served  elsewhere [25].  These  findings  support  the  change  of  the 
vaccine  strain  in  2009  from  A/Solomon  Island/3/2006  to  A/ 
Brisbane/59/2007  and  show  that,  despite  minimal  influenza 
vaccine  use  in  Kenya,  both  vaccine  strain  choices  were  appropri¬ 
ate  for  Kenya.  This  underscores  the  need  to  promote  the  south¬ 
ern  hemisphere  influenza  vaccine  uptake  in  Kenya. 

In  conclusion,  molecular  analyses  of  circulating  A(H1N1)  in 
Kenya  during  2007-2008  indicate  that  the  viruses  evolved  by 
point  mutational  drift  from  the  vaccine  strains.  These  findings 
are  similar  to  those  from  our  earlier  study  of  influenza  A  virus 
subtype  H3N2  during  the  2006-2007  season  [14].  Discovery  of 
G189N  and  parallel  mutations  Q108H  and  S125N  unique  to 
Kenya  among  A(H1N1)  isolated  in  2007-2008  highlights  the 
importance  of  continuing  surveillance  and  characterization  of 
emerging  influenza  virus  drift  variants  in  Africa. 
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